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1  Introduction 


In  this  paper  we  consider  the  monotone  linear  complementarity  problems 

Find  (x,y)  £  Rn  x  Rn  such  that  y  =  Mx  +  q,  {x,y)  >  0,  xTy  =  0,  (1.1) 

where  q  €  R”,  and  M  £  RnX".  Many  optimization  problems  can  be  formulated  as  monotone  linear 
complementarity  problems  including  the  linear  programming  problem.  We  write  (1.1)  as 

„  ,  /  Mx  —  y  +  q  \  . 

F(I'9)S(  XYe  J=°>  (L2) 

and 

(x,y)>  0  (1.3) 

where  X  —  diag(x)  and  Y  =  diag(y).  Let  S  denotes  the  solution  set  of  problem  (1.1).  A  solution 
pair  (x,  y)  is  said  to  satisfy  strict  complementarity  if  in  addition  to  complementarity  XYe  =  0,  it 
satisfies  x  +  y  >  0.  For  problem  (1.1),  we  will  be  interested  in  the  following  three  subsets  of  the 
index  set  {1, . . . ,  n) 

B  =  {i  :  x*  >  0  for  at  least  one  (a;*,  y*)  £  5}. 

Af  =  {i  :  y*  >  0  for  at  least  one  (x*,y*)  £  5}. 

J  —  {i  :  x*  =  y*  =  0  for  all  (x*,y*)  £  5}. 

It  is  clear  that  B,  A/",  and  J  form  a  partition  of  {1, . . . ,  n}. 

In  general,  for  monotone  linear  complementarity  problems  solutions  that  satisfy  strict  comple¬ 
mentarity  may  not  exist.  However,  for  linear  programming  problems,  the  existence  of  solutions  that 
satisfy  strict  complementarity  was  proved  by  Goldman  and  Tucker  [6].  In  this  case,  all  solutions 
in  the  relative  interior  of  the  solution  set  satisfy  strict  complementarity,  see  for  example  El-Bakry, 
Tapia  and  Zhang  [2]  and  we  have  J  —  0,  B  U  N  —  (1, . . . , «},  and  B  fl  Af  =  0. 

For  linear  programming  many  indicators  function  for  identifying  members  of  B  and  J\f  have  been 
proposed.  This  information  can  be  used  to  computational  advantages  in  interior-point  methods. 
A  thorough  study  for  this  issue  can  be  found  in  El-Bakry,  Tapia,  and  Zahng  [2],  In  this  paper  we 
are  concerned  with  identifying  members  of  B,  Af,  and  J .  We  consider  three  indicators  for  that 
purpose.  The  variables  used  as  indicators,  the  primal-dual  indicator,  and  the  Tapia  indicators. 
The  variables  used  as  indicators  are  perha,ps  the  most  widely  used  indicators  in  computational 
optimization.  The  primal-dual  indicators  are  also  used  for  the  purpose  of  identifying  the  above 
subgroups  of  variables,  particularly  in  the  context  of  primal-dula  interior-point  methods,  see  for 
example  Gay  [4],  Ye  [12],  and  Lustig  [8].  Monteiro  and  Wright  [10]  proposed  the  use  of  the 
primal-dual  indicator  in  the  context  of  degenerate  monotone  complementarity  problems.  The  Tapia 
indicators  were  proposed  by  Tapia  [11]  for  the  identification  of  active  and  inactive  constraints 
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for  general  constrained  optimization  problems.  The  Tapia  indicators  were  used  by  El-Bakry  [1] 
and  Mehrotra  [9].  Their  convergence  properties,  in  primal-dual  interior-point  methods  for  linear 
programming,  was  studied  in  El-Bakry,  Tapia,  and  Zhang  [2].  The  convergence  of  the  Tapia 
indicators  in  the  absence  of  strict  complementarity  has  not  been  studied. 

The  motivation  for  studying  the  Tapia  indicators  in  this  context  is  that  solutions  that  satisfy 
strict  complementarity  may  not  exist  in  general  complementarity  or  optimization  problems.  One 
example  is  the  monotone  linear  complementarity  problem  (1.1).  Another  interesting  example  is  a 
class  of  problems  arising  from  the  control  of  some  elasto-dynamic  systems,  see  for  example  [5]. 

In  this  paper,  we  prove  that  the  Tapia  indicators  corresponding  to  variables  with  indices  in  J7, 

i.e.  variables  that  do  not  satisfy  strict  complementarity  at  every  solution  ( x*,y *)  £  S,  actually 
converge.  This  result  is  a  direct  consequence  of  the  work  of  Monteiro  and  Wright  [10].  It  should 
be  noted  that  this  result  holds  only  for  feasible-iterate  algorithms,  whereas  the  convergence  of  the 
Tapia  indicators  when  strict  complementarity  holds  does  not  require  feasibility.  We  also  establish 
a  convergence  rate  of  the  Tapia  indicators  in  this  context.  Finally  we  study  the  convergence  and 
convergence  rate  for  the  variables  used  as  indicators  and  the  primal-dual  indicators. 

This  paper  is  organized  as  follows;  in  Section  2,  a  generic  algorithm  is  introduced  and  the  Tapia 
indicators  and  their  convergence  properties  for  linear  programming  problems  are  stated.  We  also 
define  two  widely  used  indicators  in  interior-point  methods,  namely  the  variables  used  as  indicators 
and  the  primal-dual  indicator.  In  Section  3  we  establish  a  convergence  result  for  the  Tapia  indicators 
corresponding  to  variables  with  indices  in  J .  We  also  study  the  convergence  and  the  convergence 
rate  for  the  three  indicators.  A  result  concerning  the  performance  of  interior-point  methods  in  the 
existence  of  solutions  that  do  not  satisfy  strict  complementarity  is  also  presented.  We  present  some 
numerical  experiments  in  Section  4.  Final  conclusions  are  given  in  Section  5. 

2  The  Algorithm  and  Indicators 

A  generic  primal-dual  Newton  interior-point  method  for  problem  1.1  has  the  following  structure. 

Algorithm  2.1  (Primal-Dual  Interior-Point  Method) 

Given  a  strictly  feasible  starting  point  (x°,y°).  For  k  =  0, 1, . . .,  do 

1.  Choose  crk  £  (0, 1)  and  set  p(x,  y)  =  rrk  ;/  . 

2.  Solve  the  following  system  for  ( Axk ,  Ayk): 

F\x\  yk)(Ax,  Ay)  =  -F(xk,yk)  +  p{xk ,yk)e  (2.1) 

3.  Choose  rk  £  (0,1)  and  set  the  steplength  ak  —  min(l,  Tkak)  where 

„k  _  _ -1 _ 

min((Afc)-1  A:rfc,  (Yk)~1  Ayk ,  —  j) 
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4-  Form  the  new  iterate 

(xk+\yk+1)  =  (xk,  yk)  +  ak(  Axk,  A  yk). 

Observe  that  since  ( x° ,  y°)  is  strictly  feasible,  all  the  subsequent  iterates  are  also  strictly  feasible. 
We  make  the  following  assumptions. 

(Al)  ( xk)Tyk  0. 

(A2)  xkyk  >  j(xk)Tyk  for  all  i  =  1, . . . ,  n  and  for  some  positive  constant  7. 

Now,  we  state  three  indicators,  the  variables  used  as  indicators,  the  primal-dual  indicator  and, 
the  Tapia  indicators  for  problem  (1.1). 

The  Variables: 

V(xk)  =  xk ,  and  V(yk)  =  yk . 


The  Primal-dual  indicator: 


where  Yk  =  diag(yk). 


PD(xk,yk)  =  ( Yk)~1xk , 


The  Tapia  indicators: 


and 


T{xk)  = 


T(yh  = 


xk+1 


■tL 

Vi 


where  xk+1  =  xk  +  akAxk  and  yk+1  =  yk  +  akAyk. 

For  the  sake  of  completeness  we  include  El-Bakry,  Tapia,  and  Zhang  [2]  convergence  result  for 
the  Tapia  indicators  for  linear  programming  problems. 


Proposition  2.1  (El-Bakry-Tapia-Zhang)  Let  the  sequence  of  iterates  {(xk,  yk)}  he  generated 
by  Algorithm  2.1  applied  to  the  linear  programming  problem.  Under  the  conditions 


1.  S  is  bounded, 

2.  ( xk)Tyk  -+  0, 

3.  6 ^  ^  7  for  k  and  some  positive  constant  7, 
4-  <Jk  — *  0  and  rk  —>  1 
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then  for  i  =  1, . . n 


0 

1 


if  i  €  M 
ifieB 


lim  T(xk)  — 


and 


Urn  T(,jf) 


0  if  i  £  B 
1  if  i  £  M 


where  xk+1  =  xk  +  (3kAx  and  yk+1  —  yk  +  /3k Ay  for  any  j3k  S  [afc,l],  where  ak  is  given  by 
Algorithm  2.1. 


An  example  that  demonstrates  the  behavior  of  the  primal  Tapia  indicator  for  AFIRO,  one  of  the 
Netlib  [3]  set  of  test  problems,  is  given  in  Figure  1. 


Figure  1:  The  Tapia  primal  indicator  for  all  variables  in  AFIRO. 


3  Convergence  in  the  Absence  of  Strict  Complementarity 

One  interesting  property  of  the  Tapia  indicators  is  that  their  0-1  separation  property  holds  component¬ 
wise,  i.e.  if  a  pair  ( xj,yj ),  j  6  M  satisfies  strict  complementarity,  then  the  0-1  separation  property 
holds  for  that  pair  even  if  other  pairs  do  not  satisfy  strict  complementarity. 

In  a  recent  paper,  Monteiro  and  Wright  [10]  considered  interior-point  methods  for  degenerate 
monotone  complementarity  problems  (problems  of  the  form  (1.1)  where  the  solution  set  <5  may  not 
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contain  a  solution  that  satisfies  strict  complementarity).  It  is  satisfying  that  as  a  direct  consequence 
of  their  result,  we  are  able  to  prove  convergence  of  the  Tapia,  indicators  for  pairs  (zj,?/j),  j  £  J . 

Lemma  3.1  Let  {(xk,  yk)}  be  a  sequence  generated  by  Algorithm  2.1  and  that  Assumptions  (Al)- 
A  (2)  hold.  Assume  further  that 

1.  the  centering  parameter  is  chosen  such  that  ak  — ►  0. 

2.  lim  ak  — *•  a  £  (0, 1]. 


Then 


Moreover 


and 


{1  —  a  if  i  £  J\f 
1  -  f  if  i  e  J 
1  if  ieB 

| T(x?)  -  1|  <  C(^(xk)Tyk  +  ak )  for  i  £  B , 
\T(xk)\  <  C(y/(xk)Tyk  +  ak)  for  i  £  J\T. 


(3.2) 

(3.3) 


Proof:  Let  i  £  fi,  then 


,Ax 


xk+1 

lim  =  lim  (1  +  ak  k 

oo  X-  t-too  X 


But  from  Lemma  2.2  and  Lemma  4.4.  of  Monteiro  and  Wright  [10], 


**  >  Cl  and  Axk  =  0([(xk)T yk}^2  +  ak)  for  i  £  B , 


where  C\  is  some  positive  constant.  Thus 


since  ( xk)T yk  — ►  0  and  ak 
tarity 

we  have 

Therefore 


lim 


-  1, 


k-+  oo  X- 

l 

0  by  the  assumptions  of  the  lemma.  From  the  linearized  complemen- 


Axiyi  +  A  yiXi  =  -x.iyi  +  axTy, 


„fc+i  „,fc+i 


,  Vi  rk,  ,  _k_k(x  )  v 

— j —  +  — J-  =  (2  -  a  )  +  a  (7  — ■ 


lim  7 .  —  1  —  a. 
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From  symmetry  we  have,  for  i  G  Af, 


_*+i 


lim 


1  —  a. 


k—foo  x “ 

From  Lemma  4.5  in  Monteiro  and  Wright  [10]  we  have 


Ax ; 


1 


Hence,  for  i  £  J,  we  have 


lim  — p-  = - for  i  G  J . 

k—too  x-  2 


+  1  .  A-.*: 

lim^oo  -jpr  =  lim^oo(l  +  «  ^) 

=  i_2 

2  ‘ 


Now  from  Lemma  4.4.  of  Monteiro  and  Wright  [10]  we  have,  for  i  6  B, 

I  A**  I 


|T,(*fc)  -  1|  <  oA 


x ^ 


<  C(sJ(xk)TZk  +  Ofc). 


(3.4) 


It  is  easy  to  see  that  the  same  result  is  also  true  for  Tt(yk )  with  i  G  W.  Now  we  have,  for  i  G  AT, 

\Ti(xk)\  <  \Ti(yk)  -  1|  +  0(ak), 


which  implies  that 

| Ti(xk)\  <  C(\J(xk)Tzk  +  crk)  for  i  G  N. 

and  completes  the  proof. 

Lemma  3.1  states  that  the  Tapia  indicators  converge  if  ak  — ►  a.  It  is  not  clear,  however,  that 
the  sequence  {a^}  converges  in  many  interior-point  algorithms  for  problem  (1.1).  Nevertheless, 
{a^}  is  usually  constructed  to  be  bounded  in  (0, 1).  So  liminfV^oo  ak  exists.  In  this  case  a  similar 
proof  to  that  of  Lemma  3.1  gives 


lim  inf  T(xk)  =  { 

k — ►oo  | 


1  —  ex'  if  i  G  Af 
l  - Y  if  iej 
1  if  i  G  B 


where  a'  =  liminf^oo  ak .  Therefore  as  long  as  there  is  at  least  one  limit  point  a'  of  {ak}  that  is 
not  very  small,  the  sequences  {T(xk)}  accumulates  around  three  reasonably  distinguished  values 
1  —  a,  1  —  f ,  and  1  corresponding  to  variables  with  indices  in  Af,  J ,  and  B  respectively.  This  shows 
that  the  Tapia  indicator  sequence  gives  a  reasonable  identification  criterion  even  if  the  steplength 
sequence  does  not  converge.  On  the  other  hand,  a  better  convergence  property  can  be  obtained  if 
the  definition  of  the  Tapia  indicator  is  slightly  modified  changed  as  follows 


T(xk) 


xk  +  Axk 
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and 


In  this  case  we  have 


0  if  i  £  Af 

lim  T(xk)  =  <  1/2  if  i  £  J 

k—+  oo 

1  if  I  e  B 


Remarks: 

•  The  convergence  rate  results  (3.2)  and  (3.3)  also  hold  for  Ti(yk)  for  i  £  Af. 

•  Assume  that  ok  =  0[(xk)Tyk]  in  Algorithm  2.1.  From  (3.2),  (3.3),  and  Theorem  2.6  of 
Monteiro  and  Wright  [10]  we  conclude  that  the  convergence  rate  of  the  Tapia  indicators,  for 
i  £  BUAf,  is  at  least  R-linear.  This  result  is  weaker  than  the  rate  obtained  by  El-Bakry,  Tapia, 
and  Zhang  [2],  As  we  will  see  in  the  numerical  experimentation,  although  the  Tapia  indicators 
demonstrate  a  slow  local  rate  of  convergence,  they  give  a  reasonable  early  identification  of 
the  status  of  the  variables. 


The  following  corollary  is  a  direct  result  of  the  Lemma  3.1  and  concerns  with  the  convergence 
and  rate  of  convergence  of  both  the  variables  used  as  indicators  and  the  primal-dual  indicator. 

Corollary  3.1  Let  {(xk,yk)}  be  a  sequence  generated  by  Algorithm  2.1  such  that  Assumptions 
(A1)-(A2)  hold.  Assume  further  that 

1.  the  centering  parameter  is  chosen  such  that  ak  —*  0. 

2.  lim  ak  — ►  a  £  (0, 1]. 

Then  for  all  i  £  Af,  both  the  variables  as  indicator  V(x)  and  the  primal-dual  indicator  P  D(x,  y) 
converges  to  zero  Q-linearly  with  Qi-factor  1  -  a,  and  for  all  i  £  J ,  V(x)  converges  to  zero 
Q-linearly  with  Q\-factor  1  —  f 

Moreover,  if  ak  — ►  1  then  for  i  £  Af,  both  V (x)  and  PD(x,  y)  converges  to  zero  Q-superlinearly, 
and  for  all  i  £  J ,  V(x )  converges  to  zero  Q-linearly  with  Qi- factor 

Proof:  The  proof  follows  from  Lemma  3.1  and  the  fact  that  the  Tapia  indicator  T{[xk)  is  the 
Qi-factor  for  both  {xk}  and  {xk/yk}. 

It  is  worth  mentioning  that  while  the  variables  used  as  indicators  converge  for  variables  with 
indices  in  J ,  the  convergence  of  the  primal-dual  indicator  has  not  been  yet  established  for  these 
variables.  The  best  that  has  been  proved  for  those  variables  is  that  their  primal-dual  indicators 
are  bounded.  This  conclusion  follows  directly  from  Lemma  2.4  of  Monteiro  and  Wright  [10].  This 
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fact  lead  Monteiro  and  Wright  [10]  to  propose  the  use  of  the  primal-dual  indicators  only  to  predict 
members  of  AT  and  B  in  their  finite  termination  technique.  An  example  of  the  behavior  of  the 
three  indicators  is  given  in  the  next  section.  The  example  actually  suggests  that  the  primal-dual 
indicators  for  variables  in  J ,  if  they  converge,  do  not  converge  to  meaningful  limits  that  can  be 
used  early  in  the  algorithm  to  detect  the  existence  of  variables  ( Xi,yi ),  i  €  J. 

Montiero  and  Wright  [10]  proved  that  the  convergence  rate  of  {(xk)T y}  is  at  most  Q-linear  if 
They  also  note  that  the  Qi-fnctor  in  this  case  is  less  or  equal  1/4  if  ak  — ►  1.  A  final  result 
concerning  the  rate  of  convergence  of  Algorithm  2.1  follows  from  Lemma  3.1,  and  is  given  in  the 
following  proposition.  This  result  states  that  the  best  Qi  -factor  that,  can  be  obtained  in  this  case 
is  1/4. 


Proposition  3.1  Let  {(xk,yk)}  be  a  sequence  generated  by  Algorithm  2.1  such  that  Assumptions 
(A1)-(A2)  hold.  Assume  further  that 

1.  J  ±  0- 

2.  the  centering  parameter  is  chosen  such  that  ak  — »•  0. 

3.  lim^oo  ak  -*■  1. 

Then  the  the  sequence  {(xk)7 yk}  converges  to  zero  Q-linearly  with  Qi -factor  l/\. 


Proof:  From  Theorem  4.1  in  Monteiro  and  Wright  [10],  the  -factor  of  {(xk)T yk}  is  given  by 

(3.5) 


Q  i  =  lim  sup 

k— i-oo 


which  implies,  under  Condition  3,  that 
From  Lemma  3.1,  we  have 


X  _  k  ,  (a*)2  Q >'kj)Tykj 
+  4  (*W’ 


Q 1  -  4‘ 


xk+1yk+1  1 

lim  — — t— r —  =  -  for  all  i  6  J . 

„,K  ^ 


k-koo  x*y* 

Consider  the  Qi-factor  of  the  sequence  {dj}  where  Oj  —  YlizJ  xiVi- 


,  V  xk+1vk+1 

Qi({9j})  =  lim  supfc^oo^W-^A 

,  2-,,£JXiyi 


(3.6) 

(3.7) 


=  lim  snpk^oo 
>  lim/t^oo  min 


Tnej  (E rtiejXj-yj) 
^  xfyf  J  4’ 


-1  **+v+* 


xiVi 


Eiej 


(3.8) 
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since  all  the  sequences  {X|  ky'k — },  for  i  6  J ,  converges  to  1/4  by  (3.7).  It  is  also  clear  from 

xi  Vi 

Lemma  3.1  that 

&BuA r  0  q- super  linearly  , 


where 

E  *.-»*• 

ieBuAf 

Now  consider  the  Qj-factor  of  {{xk)Tyk}, 


Qi({(xk)Tyk})  = 


> 


limsup^^ 
lim  supfc^^ 
lira  sup^^ 


E« 


k+ 1  k+ 1 
nxt  yi 


E« 


na4i/f 
4+i 


gy-1  ' 

9j 


/  i  + 


(3.9) 


Since  {6*^-}  converges  to 
have 


zero  Q-linearly  with  Qi -factor  at  least  1/4,  then  for  sufficiently  large  k  we 


4 


>  r  >  0. 


which  implies  that 


Therefore 


ok 

a8uAf  a 
%  ^ 


Oi({(*W»  > 


i 

4' 


From  (3.6)  and  the  last  relation,  we  obtain  the  desired  result. 

This  proposition  asserts  that,  for  degenerate  monotone  linear  complementarity  problems,  the 
best  QrTactor  of  the  sequence  {||F(a;fc, //fc))|1},  where  {(cfc,  yk)}  is  generated  by  Algorithm  2.1.  is 
1/4.  The  existence  of  even  one  pair  ( x,j,y3 ),  j  6  J  adversely  slows  down  the  algorithm  regardless 
of  how  fast  other  variables  approach  zero.  This  fact,  will  be  numerically  demonstrated  in  the  next 
section. 


4  Numerical  Example 

In  the  following  experiment,  a  perturbed  Newton  interior-point  algorithm  is  used.  The  algorithm 
is  coded  in  MATLAB.  A  problem  is  said  to  be  solved  to  an  accuracy  of  10-rf,  for  some  positive 
integer  d,  if  the  algorithm  stops  with  ||F(a:fc,  j/fc)||2  <  I0-ri.  The  problem  that  is  considered  here 
is  solved  to  an  accuracy  of  10-8.  The  parameter  ak  in  Algorithm  2.1  is  chosen  such  that  ak  —*  0. 
The  iterates  generated  by  the  code  are  not  necessarily  feasible  with  respect  to  the  linear  equations 
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Mx  -  y  +  q.  The  experiments  were  performed  on  a  Sun  SPARCstation  10  Model  41  with  128 
Megabytes  of  memory  running  SunOS  4.1.3. 

An  example  of  a  problem  with  a  solution  that  does  not  satisfy  strict  complementarity  is  the 
following 


minimize  (aq  +  Sx2  +  x3)2  +  4(aq  +  x2)2 
subject  to  1  —  aq  —  x2  —  £3  =  0 

6x2  +  4^3  -  z  j  -  X4  -  3  =  0 


(4.10) 


(xi,x2,x3,x4)  >  0 


This  is  problem  number  32  in  the  Hock  and  Schitkwoski  [7]  set  of  test  problems.  Figure  2.1 


Figure  2:  The  Tapia  indicator  for  all  variables  in  problem  #32. 

demonstrates  the  behavior  of  the  Tapia  indicators  for  the  all  variables  in  the  problem.  We  note 
that  £4  is  a  slack  variable  added  to  the  original  inequality  constraints.  The  algorithm  gives  the  point 
(6.26  X  10~5, 0,  0.999, 0.997)  as  an  approximate  solution  to  the  problems.  The  Lagrange  multipliers 
corresponding  to  the  inequality  constraints  are  (5.01  X  10— 4 , 3.999,  0, 0).  Hence  both  aq  and  the 
corresponding  multiplier  are  zero  at  the  solution.  We  note  that  for  this  example  we  have  ak  — >  1. 
The  behavior  of  the  Tapia  indicator  predicted  by  Lemma  3.1  is  clearly  demonstrated  where  the 
Tapia  indicator  corresponding  to  x\  converges  to  0.5. 

A  comparison  between  the  three  indicators  is  given  in  Table  1  for  this  problem.  The  first 
column  gives  the  iteration  number.  The  second  column  gives  the  Tapia  indicator  for  the  variable 
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Iteration  number 

Ti(*k) 

Ti(yk) 

Tk 

AhA 

1 

1.4096e+00 

1.8187e+00 

1.0000e-01 

1.9791e+01 

5.0528e-03 

2 

1.0678e+00 

6.3976e-02 

1.4096e-01 

3.5994e+01 

3.9161e-03 

3 

6.5954e-01 

5.1126e-01 

1.4882e-01 

8.2617e+00 

1.8013e-02 

4 

8.3130e-01 

3.1037e-01 

9.8153e-02 

4.2239e+00 

2.3238e-02 

5 

7.3057e-01 

3.7558e-01 

8.1708e-02 

1.3310e+00 

6.1391e-02 

6 

5.2704e-01 

5.0405e-01 

5.9693e-02 

4.9987e-01 

1.1942e-01 

7 

5.0107e-01 

5.0053e-01 

3.1461e-02 

2.5196e-01 

1.2486e-01 

8 

5.0006e-01 

5.0004e-01 

1.6002e-02 

1.2803e-01 

1.2499e-01 

9 

5.0001e-01 

5.0000e-01 

8.0022e-03 

6.4018e-02 

1.2500e-01 

10 

5.0000e-01 

5.0000e-01 

4.0084e-03 

3.2067e-02 

1.2500e-01 

11 

5.0000e-01 

5.0000e-01 

2.0042e-03 

1.6034e-02 

1.2500e-01 

12 

5.0000e-01 

5.0000e-01 

1.0021e-03 

8.0168e-03 

1.2500e-01 

13 

5.0000e-01 

5.0000e-01 

5.0105e-04 

4.0084e-03 

1.2500e-01 

14 

5.0000e-01 

5.0000e-01 

2.5053e-04 

2.0042e-03 

1.2500e-01 

15 

5.0000e-01 

5.0000e-01 

1.2526e-04 

1.0021e-03 

1.2500e-01 

Table  1:  Several  indicators  in  the  absence  of  strict  complementarity 


x\  while  the  third  gives  the  Tapia  indicator  for  the  corresponding  multiplier  y\.  The  fourth  and 
fifth  columns  give  the  value  of  the  variables  and  \j\  respectively,  while  the  sixth  column  gives  the 
value  of  the  primal-dual  indicator  xi/y^.  From  the  table  we  can  see  that  although  the  variables  X] 
and  yi  both  converges  to  zero,  none  of  them  is  “sufficiently”  small  until  iteration  10.  In  fact  both 
Xi  and  y\  are  not  small  enough  even  when  the  algorithm  terminates.  On  the  other  hand,  although 
the  primal-dual  indicator  x/y  seems  to  converge  to  0.125,  its  terminal  value  is  meaningless.  On 
the  contrary,  both  Tapia  indicators  converge,  as  the  theory  predicts,  to  0.5.  Observe  also  that  they 
converge  relatively  early  (starting  from  iteration  6).  We  note  here  that  at  iteration  6,  the  residual 
|| Mx  —  y  +  q ||  is  10-4,  i.e.  the  iterates  are  close  to  the  feasibility  region. 

Observe  that,  in  Figure  2,  the  Tapia  indicator  corresponding  to  the  zero  variable  x2  is  small  but 
not  close  to  zero.  This  is  due  to  the  fact  that  the  Tapia  indicator  corresponding  to  this  variable 
converges  to  zero  slowly,  as  predicted  by  the  theory.  Observe,  however,  that  it  gives  an  early 
identification  of  the  states  of  x2  at  the  solution.  In  Table  2  the  Tapia  indicator  corresponding  to  x2 
is  given.  The  first  column  gives  the  iteration  number  (the  code  terminated  successfully  at  iteration 
15)  for  the  last  6  iterations.  The  second  column  gives  the  Tapia,  indicator  for  the  zero  variable  when 
a  =  0(\\F(xk ,  yk)\\2) .  It  is  clear  that  the  convergence  of  the  Tapia,  indicator  to  zero  is  slow.  The 
third  column  gives  the  values  of  the  same  indicator  when  a  =  0(\\F(xk,  2/^)||4)-  The  convergence  to 
zero  is  obviously  faster.  In  fact  it  seems  that  the  Tapia  indicator  T2(xk)  convergence  to  zero  linearly 
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with  a  Qi-factor  close  to  0.50.  This  means  that  zero  components,  either  in  x  or  y,  converges  faster 
to  zero.  It  is  interesting,  however,  to  observe  that,  the  algorithm  itself  was  still  slow  (it  took  the 
same  number  of  iteration  in  both  cases).  Another  interesting  point  to  observe  is  that  in  both  cases 
the  Qi-factor  of  the  residual  sequence  is  approximately  1/4,  see  columns  4  and  7.  Other  runs  with 
tighter  stopping  criterion  also  support  our  observation.  In  all  cases  the  best  Q\ -factor  we  obtained 
was  1/4.  In  all  cases  ak  — ►  1.  This  performance  is  predicted  by  Proposition  3.1. 


Iteration 

a 

*  =  ||ir0n*,yfc)||l 

(7 

=  ||  F(x\yk 

m 

number 

relative 

Tapia 

Q 1-factor 

relative 

Tapia 

Ql-factor 

residual 

indicator 

of  residual 

residual 

Indicator 

of  residual 

3.21e-05 

1.9410e-02 

2.4884e-01 

4.90e-05 

4.9492e-03 

2.5521e-01 

8.03e-06 

1.6790e-02 

2.5016e-01 

1.24e-05 

2.48()9e-03 

2.5306e-01 

2.01e-06 

1.6458e-02 

2.5031e-01 

3.11e-06 

1.2465e-03 

2.5081e-01 

5.02e-07 

1.6060e-02 

2.4975e-01 

7.78e-07 

6.2397e-04 

2.5016e-01 

14 

1.26e-07 

1.5854e-02 

2.5100e-01 

1.95e-07 

3.1208e-04 

2.5064e-01 

15 

3.14e-08 

1.5742e-02 

2.4921e-01 

4.87e-08 

1.5605e-04 

2.4974e-01 

Table  2:  Tapia  indicator  for  two  choices  of  crk 


5  Conclusion 

The  convergence  and  the  convergence  rate  of  the  Tapia  indicators  for  feasible-iterate  interior-point 
methods  for  monotone  linear  complementarity  problems  is  established  when  strict  complemen¬ 
tarity  does  not  hold.  It  is  shown  that  the  Tapia  indicators  converge  at  least  R, -linearly  to  their 
terminal  value.  The  convergence  rate  for  the  variables  used  as  indicators  and  the  primal-dual  in¬ 
dicator  corresponding  to  zero  variables  is  also  established.  Another  result  concerning  the  behavior 
of  interior-point  methods  in  the  presence  of  solution  that  do  not  satisfy  strict  complementarity 
follows  form  Proposition  3.1.  This  proposition  asserts  that  the  best  Qi -factor  of  the  sequence 
{( xk)Tyk }  is  1/4.  The  convergence  of  the  Tapia  indicators  for  infeasible-iterate  methods  is  still 
under  investigation. 
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